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Abstract 



We formulate and study a class of U(N) -invariant quantum mechanical models 
of large normal matrices with arbitrary rotation-invariant matrix potentials. We con- 
centrate on the U(N) singlet sector of these models. In the particular case of quadratic 
matrix potential, the singlet sector can be mapped by a similarity transformation onto 
the two-dimensional Calogero-Marchioro-Sutherland model at specific couplings. For 
this quadratic case we were able to solve the TV— body Schrodinger equation and obtain 
infinite sets of singlet eigenstates of the matrix model with given total angular momen- 
tum. Our main object in this paper is to study the singlet sector in the collective field 
formalism, in the large-iV limit. We obtain in this framework the ground state eigen- 
value distribution and ground state energy for an arbitrary potential, and outline briefly 
the way to compute bona-fide quantum phase transitions in this class of models. As 
explicit examples, we analyze the models with quadratic and quartic potentials. In the 
quartic case, we also touch upon the disk-annulus quantum phase transition. In order 
to make our presentation self-contained, we also discuss, in a manner which is somewhat 
complementary to standard expositions, the theory of point canonical transformations in 
quantum mechanics for systems whose configuration space is endowed with non-euclidean 
metric, which is the basis for constructing the collective field theory. 



* e-mail: joshua@physics. technion. ac. il 
**permanent address 



1 Introduction 



Random normal matrix models were used recently in studying Laplacian growth 
processes, and in particular, the fingering instability of the boundary of various two- 
dimensional fluids [1], the theory of the r-function for analytic curves [2] and two- 
dimensional string theory [3]. 

An earlier important use of normal matrix models, closely related to the works 
mentioned above, was in quantum Hall physics [4], where it was shown that the 
partition function of the normal matrix model coincided with the zero-temperature 
partition function of two-dimensional electrons in strong (uniform, as well as varying) 
magnetic fields, and its identification as a complexified form of the Toda lattice r- 
function was established. (The electron quantum Hall liquid was one of the systems 
studied in [1].) 

The structure of correlation functions in the normal matrix models was studied 
in detail in [5, 6] and later in [7]. 

Normal matrices were also mentioned recently in [8], as a specific kind of complex 
random matrices which evaded a certain generic geometrical constraint on the shape 
of the two-dimensional eigenvalue distribution of a large class of circularly symmet- 
ric probability ensembles of complex random matrices, known as the "Single Ring 
Theorem" [9]. 

All the studies mentioned above employed time independent normal matrices. 
An interesting quantum mechanical model of truly dynamical normal matrices, with 
quadratic matrix potential, was introduced in [10], in their study of the physics of 
two-dimensional long-range Bose liquids. It was pointed out in [10] that the U (N) sin- 
glet sector of their normal matrix model was equivalent to a certain two-dimensional 
generalization of the Calogero-Sutherland model [11, 12], which contained long-range 
three-body interactions, in addition to the repulsive two-body interactions familiar 
from the one-dimensional model, at specific values of the couplings of those inter- 
actions. (For more details concerning this equivalence see section 2.2.) Under this 
equivalence, the N complex eigenvalues are mapped onto the positions of the interact- 
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ing particles of the Calogero-Sutherland model in the plane, in a manner analogous to 
the relation between the Dyson matrix ensembles and the one-dimensional Calogero- 
Sutherland model [12]. (See also section 3 of [13].) 

It should be noted that a three-dimensional version of the Calogero model, very 
similar to the one studied in [10], was formulated and partly solved long ago in [14]. 
Thus, we shall refer to these higher-dimensional versions of the Calogero model (with 
the Sutherland modification of the quadratic piece of the potential) more appro- 
priately as Calogero-Marchioro-Sutherland (CMS) models [15, 16]. We should also 
mention the construction of multidimensional Calogero models and their relation to 
matrix models in [17], as well as the generalized multidimensional Calogero models 
discussed recently in [18]. 

In this paper we shall generalize the quantum mechanical normal matrix model 
of [10] into a large class of models with U(N)- and rotation-invariant potentials. 
We will then analyze the singlet sector of these models in the large- N limit, in the 
framework of quantum collective field theory [19, 20]. We will focus on studying 
the ground state of our generic quantum mechanical model. Thus, given an arbitrary 
U (N)- and rotation-invariant potential, we will determine the ground state eigenvalue 
distribution and ground state energy explicitly, in the large- N limit. In a separate 
publication [21], we will construct and analyze the effective hamiltonian of small 
fluctuations around the ground state configuration. 

This paper is organized as follows. In the next section we discuss the quantum 
mechanics of the normal matrix model, concentrating on the singlet sector. In par- 
ticular, we solve the N— body Schrodinger equation with quadratic matrix potential 
and present infinite sets of exact singlet eigenstates with given angular momentum. 
We also discuss in this section the mapping of the singlet sector of the normal matrix 
model onto the two-dimensional CMS model. In Section 3 we will construct the col- 
lective field quantum hamiltonian corresponding to the singlet sector hamiltonian H s 
in (2.15). The standard construction of collective field hamiltonians, as presented in 
[19, 20], is formulated for systems defined in a flat configuration space, endowed with 
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cartesian coordinates. Such coordinates are not explicitly available for our system of 
normal matrices, as will be explained in that section. In general, the transformation 
to the collective hamiltonian is achieved by a quantum mechanical point canonical 
transformation from the original dynamical variables to the collective coordinates. 
Thus, we begin section 3 with a general brief discussion of quantum mechanical point 
canonical transformations in a configuration space endowed with a non-euclidean met- 
ric. We believe our presentation and results in this part of section 3 supplement those 
presented in chapter 6 of [20]. We then use these results to construct the collective 
field theory of our quantum mechanical model of normal matrices. We show that the 
Heisenberg equations of motion (as well as the classical Hamilton equations of mo- 
tion) of this model can be interpreted as the equations of motion of an eulerian fluid, 
similarly to an analogous interpretation of the equations of motion of the collective 
field theory of hermitean matrices [25] , or its equivalent formulation as fermionic field 
theory [26]. 

Finally, in section 4 we find the ground state eigenvalue distribution and ground 
state energy of this collective hamiltonian (and thus, those of the original model) 
for an arbitrary matrix potential (2.14). We also outline briefly the way to compute 
bona-fide quantum phase transitions in this class of models. As explicit examples, 
we analyze the models with quadratic and quartic potentials. In the quartic case, we 
also touch upon the disk-annulus quantum phase transition. 
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2 The Quantum Mechanical Normal Matrix Model 



Let us first recall some basic facts about normal matrices. Consider the N x N 
complex normal matrix M. Thus, 

[M,M*] = 0, (2.1) 

which means M and M' are diagonalized by the same unitary matrix U. Conse- 
quently, any normal matrix M can be decomposed as 

M = U j ZU (2.2) 

where 

Z = diag (z u ...,z N ) (2.3) 

are the complex eigenvalues. Note that U in (2.2) is not unique, since there is freedom 
to multiply it on the left by a diagonal unitary matrix of N arbitrary phases. Thus, 
the count of real independent degrees of freedom in M is the 2N real and imaginary 
parts of the eigenvalues, plus the iV 2 real independent parameters of U, less the N 
arbitrary phases, adding up to iV(iV + 1) real independent parameters. 

The geometry of the manifold of normal matrices is defined by its embedding in the 
larger space of complex matrices, with the embedding effected by the constraint (2.1). 
Thus, its line element is inherited from the euclidean line element ds 2 = Tr dM^dM 
of the latter, by substituting (2.2) for M. Therefore, writing 

dM = U\dZ+ [dR,Z])U, (2.4) 

where 

dR = UdU^ , dR ] = -dR (2.5) 

is the right-invariant form (i.e., invariant under U — > UV with a fixed unitary V), we 
obtain the normal matrix line element as 

ds 2 = ]T \d Zi \ 2 + dR*j dRij \zi - Zj\ 2 . (2.6) 

i i<j 
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Note that this line element is independent of the diagonal elements dRa (which are 
essentially the differentials of the N arbitrary phases mentioned above). It depends 
only on the differentials of the real and imaginary parts of the z^, and on the differ- 
entials of the real and imaginary parts of the elements in the lower-triangular part 
of dR, N(N + 1) in all. Note also that there are no metric elements in (2.6) which 
mix dz's and dR's: motions tangential to eigenvalues and to the unitary group are 
orthogonal. 

The laplacian associated with the metric (2.6) is thus the sum of two terms 

V 2 = V 2 + V 2 U{N) , (2.7) 

where V 2 arises from pure eigenvalue variations and V 2 /^) arises from pure U(N) 
rotations. In particular, Vfj^ annihilates U(N) singlets. Thus, V 2 = V 2 in the 
U (N) singlet sector. Clearly, V 2 and V 2 ;^ are analogous, respectively, to the radial 
and to the L 2 /r 2 parts of the laplacian in atomic physics. 

The explicit expressions for V 2 and V 2 /^ can be computed from (2.6) in a 
straightforward manner. The angular part of the laplacian is 

= 211^4. (2-8) 

and the singlet part is 



i<j \ Z i Z j\ 



v 2 = —Y 



— (|A| 2 — ) + — ( |A| 2 — ' 

dxi \ dxij dyi \ dy iy 



(2.9) 



where 



i>j 

is the Vandermonde determinant of eigenvalues, and Zi = Xi + iyi. In terms of the 

complex derivatives 

d 1 d d (91(9(9 
^ = ^ = 2 ( ^-W' ^ = ^ = 2 ( &- + W (2 - U) 

we can write V 2 more elegantly as 



In this paper we study U (N)- and rotation-invariant normal matrix models defined 
by quantum hamiltonians of the general form 

H = —-V 2 + TW(M t M) . (2.13) 

Here V is a generic polynomial potential, whose couplings scale properly with N, in 
order to ensure a well-behaved large- iV limit. We will see later-on that the required 
iV-dependence is 

y^M)^-AM!Ml, (2 , 4) 

P >i iV 

(with ^-independent g p ). 

According to general principles, the ground state of (2.13) must be a U (N) singlet. 
As was stated earlier, we will find this ground state explicitly in the large- limit. 
Thus, in this paper we shall focus exclusively on the U(N) -singlet sector of (2.13), in 
which this hamiltonian is reduced to 

tf. = ~V2 + 5>(|*| 2 ), (2.15) 
z i 

which defines the dynamics of the N complex eigenvalues Zj. It acts on singlet wave 
functions Xs( z i, z \i ■ ■ ■ , z n, z* n ) which are completely symmetric under eigenvalue per- 
mutations. 

Note from (2.12) that H s is inviariant under the interchange z i: di <-> z*,d*, i.e., 
under reflection with respect to the x-axis. This means that if \s is an eigenstate of 
H s with eigenvalue E, so is x* s i 

HsXs = E Xs ^ H sX * s = Exl ■ (2-16) 

Thus, if Xs is complex (and not just a real wave function multiplied by a phase), its 
corresponding energy will be at least doubly degenerate. 

It is clear from the metric (2.6) that H s is symmetric (Xsi\H s Xs2) = (H s x s i\Xs2) 
with respect to the measure 

N 

d»({z i , Z ;}) = \A\ 2 '[[ct 2 z k (2.17) 

k=i 
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and the inner product 

(Xsi\Xa) = f d VX* s iXs2- (2.18) 

As is evident from (2.12) (or (2.9)), all interparticle interactions are lumped into the 
kinetic term of H s , through kind of a complex gauge field 

Z-a j Zi Zj 

(where the prime indicates that the sum is over all j ^ i), and its complex conjugate 
field A* = Ai ({zji}), which multiply the first order derivatives d* and di, respectively, 

= 2]T (2d*d t + A*di + Ad*) . (2.20) 

i 

2.1 Exact Eigenstates of the Matrix Model With Quadratic 
Potential 

Let us consider now (2.15) with quadratic potential 

TrV{M^M) = l-m 2 TrM^M = ^-m 2 m 2 > . (2.21) 

2 2 i 

Remarkably, infinite (albeit incomplete) sets of eigenstates of this singlet hamiltonian 
H s can be found analytically, in a similar manner 1 to [14]. 

2.1.1 Radial Eigenstates 

It is possible to compute explicitly the ground state xo of this matrix model, and also 
an infinite set of radial excitations above it. To this end, define the radial combination 

N 

£ = EN 2 > (2-22) 
i=i 

which is symmetric under permutation of eigenvalues, and look for eigenstates of H s 
of the form x = F(q). Using the identity 

£>ftbgA= (2.23) 
i=i 1 

1 In contrast with (2.21), the quadratic potential term in [14] was oc J2i<j \ z i ~ z j\ 2 ■ 
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and a similar identity for A*, we can write the Schrodinger equation H s F(g) = EF(g) 
with potential (2.21) as 

2e — +N(N+1) _ + ( E -- m * Q y =0 . (2 . 24) 

The eigenfunction solutions of (2.24) and their corresponding eigenvalues are readily 
found as 

( N(N+1) n 1 fN(N+l) \ 

Xn{Q) = Ll 2 l \m Q )e- ms , E n = I L_I I + 2 n J m , (2.25) 

(n = 0, 1, . . .), where 

(M^i!-!) " (-l) fe (^±ii + n-l)! 

is a Laguerre polynomial. 

In particular, the nodeless gaussian eigenfunction xo = e~^ rne = e~^ m ^! 2i ' 2 is 
the ground state [15], with energy 

N(N + 1) 

E = —-^ >-m. (2.27) 

The radial states of higher energy are evenly spaced, with constant gap 2m. 

From (2.17) and (2.25), we obtain that the joint probability density of the iV 
eigenvalues in the ground state is proportional to 

|Axo| 2 = |ll(^-^)|V m ^N 2 . (2.28) 

i>j 

Starting from the two-dimensional CMS model (at couplings tuned to the point cor- 
responding to the normal matrix model), it was shown in [15] (prior to [10]) that the 
ground state wave function of the CMS model was simply 

Vo CM5 =|n(^-^)|xo, (2.29) 

i>j 

which thus led to a joint probability density for the positions of the N particles, which 
was also proportional to (2.28). It was further pointed out in [15] that (2.28) was 
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also proportional to the joint probability density of the N eigenvalues in Ginibre's 
ensemble [22] of gaussian random complex matrices, which also happens to coincide 
with the joint probability density of the gaussian (time independent) normal matrix 
model. This is, of course, in complete analogy with the one-dimensional case, in 
which the joint probability density for the positions of the N particles in the ground 
state of the Calogero-Sutherland model at specific values of the repulsive interaction, 
coincides with that of the appropriate Dyson random matrix ensemble [12]. 

It was noted in [16] 2 that the eigenstates (2.25) could be also obtained using group 
theoretical considerations. As it turns out, the normal matrix model with quadratic 
potential 3 (2.21) has dynamical symmetry 0(2, 1) 

[h, D}=ih, [h, K] = 2iD , and [K, D] = -iK , (2.30) 

generated by 

1 1 N 

1=1 

D = £E + + (2-31) 

i 

The operator h is just the kinetic part of H s , K is the conformal generator, and D 
generates dilatations. Let us define the raising and lowering operators B ± as 

B ± = -(h-K^2iD) . (2.32) 
2 

These operators, together with the full hamiltonian H s = h + K generate the SU(1, 1) 
algebra 

[H 8 ,B ± ] = ±2B ± , [B-,B + ]=H S . (2.33) 

2 This was done in section 2 of [16], in the context of the two-dimensional CMS model, and not 
directly in the normal matrix model. The adaptation of the discussion in [16] to the case of our 
normal matrix model is what follows. 

3 In the following group theoretical discussion we will set m = 1 in (2.21) for convenience (but 
will reinstate it back following (2.34)). 
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Then, the states (2.25) are obtained by applying B + repeatedly to the ground state 
Xo, namely, 

Xn = (B + Txo, (2.34) 
from which it is easy to check that indeed E n = E + 2mn . 

2.1.2 Eigenstates with Nonvanishing Angular Momentum 

The eigenstates (2.25) are rotation invariant and thus do not carry angular momen- 
tum. Here we derive an infinite set of eigenstates which do carry angular momentum. 
The derivation of these states is similar to the derivation of the radial states. 

We shall look here for eigenstates of the form F(g) g(z±, . . . z^), with g(zi, . . . z^) 
a totally symmetric holomorphic function. Using the identity (2.23) (and the corre- 
sponding one for A*), we can write the Schrodinger equation as 



2 e ^ + N(N + l)^ + (E-^ B )F 



dg Y V TV ^ J 

(2.35) 

From the identity J2i = we see that if we choose the symmetric function g such 
that dig will be independent of the index i, and furthermore, if we choose it to be 
homogeneous of some degree p such that J2i(zidig) = pg, then (2.35) will imply that 

2g— + (N(N + l) + 2p) ^-+^E--m 2 g)F = 0, (2.36) 

which coincides with (2.24) up to shifting N(N + 1) by 2p. The two restrictions on g 
uniquely fix it (up to a multiplicative constant) as g(zi, . . . zn) = (J2i z i) v ■ Thus, by 
comparing (2.36) and (2.24) we see that 

f N(N+l) 

X % ) = (z 1 + ... + z N ) p li 2 +P \mg)e" me (2.37) 
is an eigenfunction of H s , with corresponding eigenvalue 

e (p) = ^(iv + 1) +p + 2n ^ m (2 3g) 

(n = 0, 1, . . .). The radial eigenstates (2.25) are just the states (2.37) with p — 0. Of 
course, x$* * s a ^ so an eigenstate, linearly independent of xffl (f° r P ^ 0), and with 
the same eigenvalue E&\ in accordance with (2.16). 
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By applying the angular momentum operator 

L z = YJL*idi - z*d*) (2.39) 

i 

to Xn^ we see that ^ carries angular momentum p. Thus, p is a nonegative inte- 
ger. Similarly, xffl* carries angular momentum —p. It is reasonable to expect that 
analogous states should exist also in the two-dimensional CMS model at arbitrary 
couplings. 

The generalized mult i- dimensional Calogero model discussed recently in [18] allows 
for N distinguishable particles in arbitrary number of dimensions. Part of the spec- 
trum of that model was studied using group theoretical methods. In two dimensions, 
and for the special case of identical particles, the model in [18] can be mapped by 
a similarity transformation onto our singlet sector hamiltonian (2.15) with quadratic 
potential (2.21). (See the discussion following (2.46).) Our eigenstates (2.37) coincide 
with a certain subset of the eigenstates obtained in [18] upon tuning the couplings of 
the latter model to coincide with ours. 

By comparing (2.38) and (2.25), we see that the pair of first excited states among 
the states we have constructed is the one corresponding to n = and p — 1, with 
energy difference 

_ = m (2.40) 

above the ground state, which thus sets an upper bound on the gap of excitations 
in the spectrum of H s . This upper bound on the gap, obtained from exact singlet 
eigenstates, is half the gap computed recently in [23] in the large-iV limit, using 
collective field analysis, for the equivalent CMS model at the normal matrix point in 
coupling space. 

2.1.3 The Quantum Normal Matrix Model and Physics of Electrons in 
the Lowest Landau Level 

The joint probability density (2.28) for N eigenvalues in the ground state xo coincides, 
upon identification of -4= as the magnetic length, with that of N noninteracting elec- 
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trons which occupy the lowest Landau Level (LLL) (at minimal angular momentum). 
Indeed, Axo can be identified with the v — 1 Laughlin wave function. 

However, this relation between the ground state of H s with quadratic potential 
and many-body wave functions in the LLL, does not persist to the excited states 
Xn (n > 1) and Xn^* (P — 1) i n (2.25) and (2.37), for a very simple reason. Our 
quantum mechanical normal matrix model is invariant under reflection with respect 
to the x-axis, which leads to the degeneracy (2.16). In contrast, in the quantum 
Hall problem, iV-electron states in the LLL (in the symmetric gauge) are in the form 
of an antisymmetric holomorphic function /a(zi, • • • , zn) ( multiplied by the gaus- 
sian e~^ me ). Any dependence in Ja on nonhomolorphic variables z* means that the 
corresponding state belongs to higher Landau levels 4 . Thus, for example, antiholo- 
morphic wave functions will correspond to higher energies than their holomorphic 
counterparts, unlike (2.16). As a particular example, consider the pairs of degenerate 
eigenstate Xo\ Xo^* °f H s . The wave function is of course a linear combination 

of LLL states 5 , but Axo*'* is of higher energy in the Landau problem. 

In very intense magnetic fields, which in our normal matrix model corresponds to 
large m 2 , the gap from the LLL to the next Landau level becomes very large, and 
low energy physics can be described entirely in terms of antisymmetric holomorphic 
wave functions. This means that the number of degrees of freedom is cut in half: to 
each electron corresponds the combination zi = Xi + iyi of position operators, and 
not Xi and yi separately. However, in our inverstigation of the singlet sector of the 
normal matrix model in this paper we do not have any special interest in taking the 
limit m — > oo so as to make the gap to the excited states infinite. Thus, we will have 
to take into account all singlet states, holomorphic as well as nonholomorphic. 

In fact, it might be the case that a version of the quantum normal matrix model 

4 As is well-known, this happens because the Landau hamiltonian in the symmetric gauge is the 
sum of a rotation invariant piece and a piece proportional to BL Z = Bj2ii z idi — z*d*), which breaks 
the symmetry under interchanging Zi,di <-» z* , d* . 

5 The p = 1 case is special, since can be written as a single determinant (rather than a 

linear combination of determinants) of holomorphic monomials z n , multiplied by gaussian factors, 
of the form Axo, where A differs from the Vandermonde determinant by having its last row entries 
zf instead z™~ x . 
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considered in this paper, in which a random piece is added to the potential \/(|,2| 2 ), 
could teach us something relevant for understanding the levitation (or floating) prob- 
lem of the extended states away from the center of the Landau band in quantum Hall 
physics, at weak magnetic fields. 

In order to relate our quantum normal matrix model and the time-independent 
normal matrix models [1, 4], however, the large m 2 limit is important, precisely 
because of the large energy gap, proportional to m, separating the ground state xo 
and the excited states of H s . Thus, at large m 2 , the low-energy physics of H s (with 
V(|,2| 2 ) = \m 2 \z\ 2 ) is captured entirely by the ground state. As was mentioned earlier, 
according to [4], the joint probability density of N noninteracting electrons, which 
occupy the LLL in a strong (and possibly nonuniform) magnetic field, coincides with 
the joint probability density of the eigenvalues of the zero-dimensional normal matrix 
model, which is given by 

P(zi, ...,z N ) = ± I - z s ) 2 e -E^(N a ) , (2.41) 

where H^(|z| 2 ) is the magnetostatic potential, namely, B(z,z*) oc 9(9*W^(|z| 2 ). This 
result was obtained by observing that the iV-electron wave function, with all spins 
parallel and with minimal angular momentum ( JV ^~ 1 - > ), was simply 

i>LLL = U( z i - h) e-^ W( ^ 2) , (2.42) 

i>j 

i.e., a generalized form of the v — 1 Laughlin wave function. For a gaussian W = 
^m 2 \z\ 2 , which corresponds to a uniform magnetic field, ipLLL is just the function 
Axq, as we saw above. 

Note, however, that for a polynomial H^(|z| 2 ) of higher degree, = e - ^^ ^N 2 ) 
is not an eigenstate of H s (2.15) with a local potential \/(|2;| 2 ). Thus, in nonuniform 
magnetic fields, we cannot relate i/jlll to the ground state of H s with any reasonable 
potential. 

As the final comment in this section, note that the many-body wave function 
Xn^ is somewhat reminescent of the single particle wave function in the nth Landau 
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level and with angular momentum p, which is given (in the symmetric gauge) by 
z p {mg)e~^ me (where in the latter expression g = \z\ 2 ). 

2.2 Relation to the Two-Dimensional Calogero-Marchioro- 
Sutherland Model 

The hamiltonian of the CMS model is [14, 15] 

h cms = -^E^ + ^E^ + GEE'TT^ + yE 1 '? (2-43) 

i i<j ij i j<k L ij L ik * i 

where = — rj. Thus, Hqms contains a three-body long-range interaction, 
in addition to the repulsive two-body long-range interaction familiar from the one- 
dimensional Calogero-Sutherland model. 

In order to see the equivalence between our normal matrix model and the two- 
dimensional CMS model (at the specific point in coupling space) we have to remove 
the "gauge fields" Ai and A* (2.19) from (2.20). These fields can be gauged away 
by a nonunitary gauge transformation, or more precisely, by a (singular) similarity 
transformation. Indeed, it is a straightforward calculation to show that 



Thus, 

I I \ i i i 



A A* 



(2.45) 



where we used Ad* ft = S^i + ttt = V?. The hamiltonian obtained from H s after the 
similarity transformation, 

H, = | A|ff * = 1 £ V ? + £ V(N 2 ) + \ E 8 -t 9J §- . (2-46) 

I I i i i 

has the standard kinetic term. Thus, we have removed the gauge fields in (2.20) at 
the price of introducing the interaction terms J2i a~~a*~' 
In order to proceed we note, using (2.19), that 
8 A3* A* 1 1 

A A * ^— -J ^—^ -v -v ^— -J 



a a* vyzi-zj-t z t - 4 
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i<j \ Z i Z j\ Z i j<k \ Zi Z j\ 2 \ Z i Z k\ 2 

= 2E4 + 2EE'^, (2-47) 

i<j ij i j<k L tj L ik 

where we denoted the position operator of the ith eigenvalue Zi = Xi + iyi in the plane 
by Yi = x{k + y{y . Substituting (2.47) in (2.46) and comparing it with (2.43), we 
see that for the particular potential | ^ | 2 ) = |m 2 |z| 2 , the hamiltonian Ti s coincides 
with H C ms at g = G = 1 [10, 15, 16]. 

The CMS model (2.43), at the special point g = G = 1, was studied recently 
in [23], in the collective field approach in the large- N limit, and it was found that 
the particles condensed into a uniform disk, in accordance with the findings of [15]. 
The authors of [23] went one step farther, and also computed fluctuations around the 
uniform ground state configuration, which required taking the first subleading term 
in the large- N expansion. They have found that the spectrum of those fluctuations 
had an energy gap equal to m, corresponding to setting V(|,2| 2 ) = \m 2 \z\ 2 in (2.46), 
as was already briefly mentioned following (2.40). In this work we extend the leading 
large- N results of [23] to the case of an arbitrary external potential. 

The hamiltonian H s acts on singlet wave-functions Xs, which are completely sym- 
metric under permutations of the Z{. The similarity transformation on H s in (2.46) 
implies a similarity transformation Xs — ► | ^\Xs( z i, z* . . . , zn, z* n ) into wave functions 
which are also completely symmetric under permutations of the eigenvalues. Thus, 
7i s in (2.46), with its conventional flat space kinetic term, corresponds to a system 
of interacting bosons in two dimensions. In contrast, in the corresponding quan- 
tum mechanical hermitean matrix models, the analog of 7i s acts on one-dimensional 
non- interacting fermions, as was shown in [24]. 

Note, however, that the similarity trasformation (2.46), which maps the matrix 
model hamiltonian H s onto the CMS hamiltonian 7i s is highly singular, since |A| = 
Uiyji^Zi — Zj) has branch points wherever two eigenvalues coincide. This might lead 
to some delicate issues in trying to decide whether the singlet sector of the normal 
matrix model and the CMS model are completely equivalent or not. There is no such 
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problem in the mapping of the singlet sector of the quantum mechanical hermitean 
matrix models onto one- dimensional noninteracting fermions. 

We stress that our construction of the quantum collective field formulation of the 
singlet sector of the normal matrix model in section 2.2 is based entirely on the non- 
euclidean configuration space with metric (2.6) and the corresponding hamiltonian 
H s . Nowhere in our construction do we resort to the singular similarity transformation 
into euclidean configuration space and CMS hamiltonian H s . Nevertheless, as we have 
already remarked above, the results of our collective field analysis (for the particular 
potential V(|z| 2 ) = \m 2 \z\ 2 ) agree (to leading order in j^) with those of [23], which 
were obtained starting with TL S . 

2.3 A Concluding Remark Concerning Cartesian Coordinates 
and Feynman Diagrams 

We end this section with the following comment: It is customary to refer to (2.2) 
as the polar decomposition of M. Thus, the Z; L and Rij (i < j) are referred to as 
the polar coordinates of M. The matrix elements are known as the cartesian 
coordinates. Thus, (2.6), and (2.8), (2.9) express the metric and the laplacian in 
polar coordinates. There are no simple expressions for these objects of normal matrix 
theory in terms of the cartesian coordinates My, since it is rather complicated to solve 
the constraint (2.1) in these coordinates explicitly. To this end one has to work in the 
bigger space of complex matrices and introduce a matrix Lagrange multiplier field 
to impose (2.1) on the dynamics, and then follow Dirac's procedure for constrained 
quantization. Interestingly enough, these complications do not arise for hermitean 
matrices, which is of course a subclass of normal matrices. In this case, the linear 
constraint M, tj = M* t implies simply V 2 = £ij d 2 / (dM^dM*^) = Tr(d 2 /dM 2 ). 

Having a lagrangian formulation of the normal matrix model in cartesian coor- 
dinates is necessary in order to formulate the large-iV "double-line" diagrammatic 
expansion of the model. In this respect, it should be noted that this expansion will 
be in powers of j? rather than powers of -r^, as should be clear from the exact ex- 
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pression (2.27) for the ground state energy in the case of a quadratic potential. 
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3 Construction of the Collective Hamiltonian 



The transformation from the hamiltonian in its original form to its expression in 
terms of collective coordinates is achieved by a quantum mechanical point canonical 
transformation from the original dynamical variables to collective coordinates. 

The configuration space of our quantum mechanical normal matrix model is 
parametrized by the polar coordinates Zi and Rij(i < j), and is endowed with the 
non-euclidean metric (2.6). 

Thus, in order to prepare the ground for the construction of the collective field 
hamiltonian for the normal matrix model, we will first recall some basic facts about 
point canonical transformations in configuration spaces whose geometry is defined by 
a non-euclidean metric, thus extending the discussion in chapter 6 of [20]. 

3.1 Point Canonical Transformations in Non-Euclidean Con- 
figuration Space 

We shall start by stating some elementary facts, mainly in order to introduce nota- 
tions. We make our points starting in Eq.(3.14). 

Consider a quantum mechanical system, defined in a D-dimensional configuration 
space C, with coordinates q a (a = 1, . . . D) and line element ds 2 = g a b(q)dq a dq b . The 
hamiltonian is given by 6 



is the invariant laplacian, expressed in terms of the coordinates q a . In (3.2) g ab is the 

inverse metric and g = det g a b. This laplacian (and the corresponding hamiltonian 

H q ) are symmetric (ipi\H q ip2) = (H q ipi\ip2) with respect to the measure y/gdq and 

the inner product (V'llV^) — / (q) ^(q)- 

c 

6 In this subsection we will write all factors of % explicitly, in order to identify the purely quantum 
mechanical terms induced by the transformation in the effective hamiltonian, defined later in this 
section. 





where 




(3.2) 
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Let us consider a point canonical transformation 

q a ^Q a = g a (q) , (3.3) 
and assume its inverse exists as well: 

q a = q a (Q) . (3.4) 

We know from rudimentary differential geometry, that since the line element ds 2 is 
invariant under this transformation, the metric in the new coordinates is given by 



dq m dq r 

0>&(Q) = ffmn(q(Q)) dQa dQb , (3.5) 



its inverse by 



and 



Thus, 



where 



^(Q)=^"(q(Q))|^^, (3.6) 
n(Q) = detn ab = g(det(^\) . (3.7) 



Vti = ^gJ, (3.8) 



is the Jacobian of the transformation, rendering the measure invariant 

VttdQ = ^/gdq. (3.10) 
Finally, the hamiltonian in the new coordinates is written 

H Q = -^V 2 Q + V(Q), (3.11) 

where 

1 ° f n ab Vn^-) (3.12) 



is the invariant laplacian in the new coordinates Q, and V"(Q) = V(q(Q)) (i.e., the 
potential is a scalar function under (3.3)). 
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Wave functions also transform as scalars, ip(q) — > ip(Q) = V ; (q(Q)) ) since ampli- 
tudes must remain invariant against the coordinate transformations: 

(iPM = J v / firfQVi*(q(Q))^ 2 (q(Q)) = / v^dq^(q)^ 2 (q). (3.13) 

c c 

The form (3.11) of the transformed hamiltonian, and the fact that wave functions 
are scalars, guarantee the invariance of matrix elements (ipi\H\ip 2 ) under coordinate 
t r ansf ormat ions . 

Let us assume that for one reason or another, the Q-coordinates represent the 
physical picture more transparently, and thus, working with Hq has some extra merit. 
Hq is symmetric with respect to the measure \^UdQ. Thus, part of the price to be 
paid working with the Q-coordinates is the need to drag along that pesky measure 
everywhere. Life would be simpler if we could rid ourselves from that measure and 
map Hq onto an effective hamiltonian H e ff which is symmetric with respect to the 
flat measure dQ. This we can achieve by performing a similarity transformation. 
Observing that 

(V>i|#|V 2 ) = / VndQri(Q)H Q MQ) = J dQ (n* Vi*(Q)) (n*H Q n-*) (q^ 2 (q 

c c 
it is obvious that required transformation is 

H eff = &H Q n- 1 * (3.14) 

with the appropriate transformation 

^(Q) = fl*^(Q) ( 3 - 15 ) 

of wave functions. Let us now massage H e ff into a more transparent form. It is a 
matter of straightforward calculation to show that 

Qab 4* - (4^^ C A ~ C « ^ ( 3 - 16 ) 



dQ a dQ b \dQ 
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where we have defined 

C a = i(logfi), a (3.17) 

and where (-) ja indicates a derivative with respect to Q a . The operator £1^ Vgf2~3 
is manifestly symmetric with respect to the flat measure dQ, as is evident in each of 
the lines in (3.16). Thus, H e ff is indeed the desired hamiltonian we set out to find, 
which, following (3.16), we may write explicitly as 

H eff = l -p a n ab p b + ^c a n ab c b + v(Q) + y (n ab c b ) a (3.18) 

where we introduced the momentum operators 

p " = -< h w- (319) 

The terms in (3.18) quadratic in h may be thought of as a generalization of the 
centrifugal barrier which arises in the radial hamiltonian in D dimensions 7 . Evidently, 
these terms are purely a quantum mechanical effect. 8 
It is easy to see that 

C a = \T b ba (3.20) 

where T^ c is the second Christoffel symbol (i.e., the connection) of il ab . However, 
sometimes a direct computation of the C a from their definition (3.17), or from the 
identity (3.20), may be too difficult to carry in practice. Thus, in order to bypass 
these potential difficulties, we shall now derive an identity satisfied by the C a , from 
which we could compute them with somewhat less effort. 

To this end we argue as follows: The invariant laplacian acting on a function which 
is a scalar under coordinate transformation produces yet another scalar function. 
Thus, 

VjlK?) = V^(?(Q)) • (3.21) 

7 The radial part of the D-dimensional laplacian = r~( D ~ 1 ' , d r (r D ~ 1 d r ), defined with respect to 
the measure r D ~ x dr, may be transformed by a similarity transformation into = rTV^r ~ = 
— — ^ = d^. — C r Cl rr C r , which is defined with respect to the flat measure dr. 

8 As a side remark, we also mention that the last term in (3.18) may be written as 4j- (fl ah Cb) = 

-afni (v 2 Q iW) . 
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In particular, the coordinate functions themselves are scalars (their differentials are 
one-forms). Let us define the quantities 

uj a = -hV 2 q Q a = -hV 2 Q Q a . (3.22) 

It follows from the definitions (3.22) and (3.17) that 

u a = -nv 2 Q Q a = (n"Vn) = -hi ab (io g fi), 6 - nn ab b 

or 

u a + 2h Q ab C b + h Q ab b = , (3.23) 

which is the desired identity to determine the C a . To simplify the computation, we are 
free to choose in (3.22) the coordinates q a in which the computation of uo a = —%V 2 Q a 
is as simple as possible. 

This concludes our brief review and exposition of point canonical coordinate trans- 
formations in quantum mechanics, which we will use in the second part of this section 
to construct the collective hamiltonian of our normal matrix model. 

The discussion of point canonical transformations in section 6.1 in [20] was presented' 
for the special case of flat euclidean metric g a b = 8 a b- For this metric, it is a straight- 
forward exercise to show that u a , C a and Q ab which we defined here for an arbitrary 
metric, coincide with their counterparts in chapter 6 of [20]. In addition, in that met- 
ric, our expression (3.18) for H e ff an d the identity (3.23) reduce to their counterparts, 
Eqs.(6.21) and (6.20), respectively, in [20]. 

3.2 The Collective Field Hamiltonian for Normal Matrices 

The ground state, as well as other £7(./V)-singlet eigenstates of the hamiltonian (2.13), 
are totally symmetric functions Xs( z i, z i, ■ ■ ■ , z n, z*) of the eigenvalues. As we saw in 
the introduction, these eigenvalues comprise a two-dimensional system of interacting 
bosons. 

9 We should mention that there is a typographical error in Eq.(6.19) in [20]. Its correct form is 
the sum of anticommutators H e f / — = j{u) a + 2h £l ab Cb + fi Q ab h , Pa}- Also note, that the last 
(total divergence) term in our (3.18) is missing from the analogous equation, Eq.(6.21), in [20]. 
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Let 



(3.24) 



be the position operator of the ith eigenvalue Zi = Xi + iyi in the plane. The density 
operators 



form a continuous basis (parametrized by r) of commuting operators which are sym- 
metric in the operators r^. Clearly, any operator which is purely a symmetric com- 
bination of i-j's may be expressed in terms of the density operators (3.25). Thus, it 
is conceivable that any singlet wave function Xs in our model could be expressed as 
functional of the eigenvalues of these density operators, namely, as a functional of 
the density function, which we also denote by <f)(r). The density function is evidently 
positive, 



This motivates us to apply the method of coordinate transformation which we 
derived in the previous subsection, and make a transformation from the eigenvalues 
Zi (or the corresponding commuting position operators r^), which live in a configu- 
ration space with the non-euclidean metric (2.6), to the set of commuting density 
operators. The latter are the collective variables for our system of bosons - they 
are the combinations of the original dynamical degrees of freedom, the r,, which are 
invariant under the permutation symmetry of the problem. The transformation we 
will perform affects only the singlet degrees of freedom, and leaves the [/(A r )-angular 
variables unchanged, of course. 

In this section we will make the transformation to these collective variables and 
study the system of bosons in the large- iV limit. This method of transformation, 
which is commonly known as the collective field method, was introduced in [19] (for 
a review, see chapter 7 of [20]). It should be thought of as a natural extension 
of the Bohm-Pines theory of high density plasma oscillations [27]. An important 
application of this method in [19] was to study the singlet sector of the quantum 



JV 1=1 




(3.25) 




(3.26) 
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hermitean matrix model. As was mentioned earlier, the same problem was solved 
earlier in [24], by mapping it onto a one dimensional gas of noninteracting fermions. 
Here, we will extend these works to normal matrices. 

The number of degrees of freedom in the singlet sector of (2.13) is 2N. In contrast, 
there is a continuum of density operators 4>(r), which are not all independent. For 
example, it is obvious that the constraint 



J dr <fi(r) 



(3.27) 



should hold. As our independent collective variables we choose the set of 2N Fourier 
modes 

r 1 N 

k = / dr e~ ikr 0(r) = - £ e"*"" (3.28) 

J iV i=l 

cutoff by k max , |k| < k max , where k is properly discretized, e.g., by putting the bosons 
in a large box of linear size L, and imposing appropriate boundary conditions. The 
details of this discretization are not important for our discussion of the large- N limit. 
Next, we assume that the eigenvalues condense into a lump of linear size L c , where 
L c « L. The average bulk density of particles in the system is thus 

N 

P- T2 , (3-29) 

and the microscopic interparticle distance will be of the order I ~ 1/yfp- Thus 
the maximal Fourier components relevant to our problem should be of the order 
kmax — 1/1 ^ yfp. We shall take the limit N — > oo together with L — > oo such 
that p » 1, or more precisely, L » L c » I. Thus, the high density limit makes 
kmax — > oo, and letting L — > oo makes k continuous. 

In this limit, which we shall assume throughout our analysis, we can consider the 
continuum of 0(r)s as independent variables, constrained by (3.27) (i.e., by 0k=o = !)■ 

To summarize, we wish to transform the normal matrix hamiltonian (2.13) from 
the coordinates q a = r^Rij and metric g ab given by (2.6), to the coordinates Q a = 
0k, Rij and metric f2 a & given by (3.5), and to construct the effective hamiltonian 
(3.18), according to the formalism developed in the previous subsection. 
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Since this transformation does not affect the U(N) -angular coordinates Rij and 
does not mix them with the eigenvalues, £l ab will have the same block-diagonal struc- 
ture as g ab , and of course, their pure £/(iV)-angular sectors will coincide. We need 
only compute the transform of the pure eigenvalue sector Q^'^' . 

We start by computing the latter, which we denote more conveniently by f2(k, k'; [(f)]), 
with the functional dependence on indicated explicitly. Since (2.6) is equivalent to 



ds 2 = J2dr 2 t +2j2 dW tJ dRij \n - rtf , 

i i<j 



(3.30) 



we obtain from (3.6) and (3.28) that 



net i / !M\ ^ k k ' k ' 
n(k,k;M)=E^-~ = 



N 



->k+k' 



(3.31) 



The effective hamiltonian (3.18) contains also the quatities C a , which we will deter- 
mine from the identity (3.23). Thus, we have to compute u^ k = u;(k; [0]). From the 
definition (3.22) we obtain 10 



w(k; [0]) 



-V 2 k 



-V 2 k , 



(3.32) 



where V 2 was defined in (2.7). Thus, from (2.9), we find 11 

1 



W (k;[0) = 



E 



-f|A| 2 ^U^f|A| 2 ^ 
dxi V / % V % 



ik 



-ikr 



. 9 log |A| 



From (2.10) we obtain 



d log |A| 



i v x » •■3) 

The latter sum has a simple physical interpretation in terms of two-dimensional 
electrostatics. If we think of the eigenvalues located at the points r k in the plane as 



10 From this point throughout the rest of the paper we set % = 1. 

11 This expression for w(k; [</>]) is the two-dimensional analog of Eq.(7.55) in [20] and Eq.(5.6) 



in [19], namely, the expression u>(k; [<fi]) = Nk 2 J da4>ku4>k(i-a) 

o 

corresponding to the hermitean matrix model with eigenvalues A, 



k 2 <f> k 



2ik 
N 



-ikX, 



i 
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positive point unit charges, then 



E^E^ (3.33) 



is the electric field at r due to all the other charges, except the one at v { . Thus, 

2ik 



W (k; [0]) = k 2 k + 1 r r -E E (i) (r.) ■ (3.34) 



We would like to avoid the i-dependence of the vector field E^(r) in (3.34) in order 
to simplify the summation. The aim is to replace it by the full electric field 



r — r,- 



However, since E(r) is singular at we have to replace E^(r) by a regulated form of 
E(r), which will tend in the large density limit to the macroscopic smoothed electric 
field, which corresponds to the smoothed macroscopic density 0(r), namely, 

E(r) = Nj dv'-^y^iv'), (3.36) 

and satisfies Gauss' law 

V • E(r) = 2vriV0(r) . (3.37) 
That (3.36) satisfies (3.37) is a trivial consequence of the identity 



r 



V • — = 2n5(r) . (3.38) 



r 



To identify the regulated form, note first that 



with s a small real parameter (which we shall assume is much smaller than the mean 
interparticle distance ~ l/^/p), should have a negligible effect on E^(r) in a typical 
configuration of the eigenvalues, in which it is most likely that |rj — rj\ > 1/y/p for 
j 7^ i, due to electrostatic repulsion. In E^(rj) we can obviously relax the restriction 
j ^ i in the summation. Thus, it is clear that 

e.M = E (^^ (3 ' 40 > 
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is the desired regulated form of the electric field, with E^(rj) = E s (rj) to be substi- 
tuted in (3.34). 

The electric field E(r) = of a point unit charge at the origin is regulated in 
this way to E s (r) = -pr^z, such that V • E(r) = 2-k5(y) — > V • E s (r) = 27r5 s (r), where 
2t[5 s {y) = ( - r 2 2 | !g2- ) 2 ■ As s tends to zero, 5 s (r) tends to 5(r). As another motivation for 
introducing the regularization (3.40), note that 

Esx ~ lEsy ~ \ \z- Zj \* + s* - Tr (,-M)t(,-M) + ^ ' (3 ' 41) 

where M is a normal matrix with eigenvalues Zi, is the regulated form of the Green's 
function G(z,z*) = Tij^ one would obtain from the method of hermitization of 
[28]. 

E s (r) is the interpolating regulated electric field which should tend, in the large 
density limit to a smooth asymptotic distribution, which we can obtain by solving 
(3.37) for E(r), with the corresponding limiting smooth distribution of <p(r) as the 
source term. Thus, in the large density limit we can replace E^(rj) in (3.34) by 
E s (rj) and write 

2?k 

W (k; [0]) = k 2 k + — • £ e-*"" E.( ri ) . 

i 

Henceforth we will suppress the index s in E s (r), with the limit s — > understood as 
the last limit taken. 

Since E(r) is of order N, as can be seen from (3.35) (or (3.36)), we should neglect 

the subleading O(N ) term k 2 k in the last equation. Therefore, to leading order in 
j_ 

N ' 

9?'k 

"(k; [fl) = E(r,), (3.42) 

i 

which is a quantity of order N. 

Thus, we have determined Q(k, k'; [0]) and u(k; [<p]). The last ingredient needed 
for determining C<^ k = C(k; [(f)]) from (3.23) is the divergence fi ab 6 . Happily enough, 
it vanishes 
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as can be seen from (3.31). Substituting (3.31), (3.42) and (3.43) in (3.23), we obtain 
the equation 

w(k; [0]) + 2£ tt(k,-k'; [</>]) C(k'; [0]) = (3.44) 

k' 

for C(k; [0]). In the combined large density and large volume limits discussed above, 
the k-sums tend to Fourier integrals. Thus, in the limit 

r H~k' 

W (k; [0]) + 2 j ^ fl(k, -k'; [0]) C(k'; [0]) = , (3.45) 

which we will transform now to r-space. 
To this end we need 

W (p; [<j,]) = J — e *- w (k; [0]) = 2V • (0(r)E(r)) , (3.46) 

where we used 

/ ^ £ e lk(r - ri) E(r,) = 1 £ 5(r - r,)E(r) = 0(r)E(r) . 

We also need 

fi(r ' r ' ; = / elk ' r+lk '' r ' fi(k ' k ' ; m = ^ Vr ' Vr ' (0(r)5(r ~ r ' )} ' (3 ' 47) 

Using (3.46) and (3.47) we Fourier transform (3.45) to r-space and obtain 

"(r; [0]) - • (0(r)VC(r; [0])) = 

|V • [0(r) (iVE(r) - VC(r; [0]) )] = . (3.48) 

The minimal solution of this equation is simply 

VC(r; [0]) = iVE(r) . (3.49) 

Thus, from (3.36), C(r; [0])) = NU(r), where U(r) is (minus) the electrostatic po- 
tential 

U(t) = X - J dr' log (r - r') 2 0(r) . (3.50) 
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VC(r; [0]) in (3.49) coincides with the leading large- N behavior of the analogous 
quantity in [23] . Also, it is the two-dimensional analog of the corresponding quantity 
obtained in the collective field formulation of the hermitean matrix model [19], 

dC{x; [0]) = N r <t>{y) dy 
dx J x — y 

which reproduces the results of [24] correctly 12 . 

In two dimensions, however, we have the freedom of adding a piece to the 
right hand side of (3.49), where u is an arbitrary divergence-free vector 

V-u = 0. (3.51) 

We shall not pursue this possibility here. 

We now have all the ingredients needed for constructing the effective collective 
hamiltonian (3.18). From (3.47) and (3.49) we obtain the collective potential V co u = 
^C a C\y as 

V coll = \J dv dv' C(r) n(r, r') C(r') = ± J dv 0(r) (VC(r)) 2 
N r 

= y /rfr0(r)E(r) 2 , (3.52) 

which is a quantity of order N 3 . In the collective field formulation of hermitean 
matrices, this term can be written equivalently as an integral over <f> 3 (x) [19, 20]. 
Similarly to (3.52), the kinetic term \P a Vl ab P b in our model is 

K=\f drdr n(r) fi(r, r') n(r') = ^- f dr Vn(r) • 0(r) Vn(r) , (3.53) 

where Il(r) is the momentum operator conjugate to 0(r). 0(r) and n(r) satisfy the 
equal-time canonical commutation relation 

[0(r),n(r')]=^(5(r-r')--^) , (3.54) 



12 In this context, we should also mention that Eq.(3.46) (which leads, together with (3.47), to 
(3.49)), is essentially the two-diemnsional analog of of Eq.(7.58) in [20] and Eq.(5.9) in [19] for 
ui(x; [</>]) in the hermitean matrix model. In those references, the one-dimensional analog of of (3.46) 
was obtained directly from the expression for uj(k; [0]), which contained the subleading k 2 4>k term. 
However, the latter term was effectively lost by taking the continuum limit of the Fourier sum into 
a Fourier integral. 
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where the inverse area subtraction arises because the zero-mode 0k=o is non-dynamical, 
due to the constraint (3.27). 

Next, we have the divergence term | (£l ab C^j = — \(uo a + n°^ b ) ta , where we used 
(3.23). Thus, from (3.46) and (3.43) we find it as 

VW- = " \ J * ^f 1 = ~\ J * V • (0(r)E(r))) . (3.55) 

It is evidently a singular boundary contribution. Assuming boundary conditions such 
that 0(r) vanishes at spatial infinity faster than r -2 , we obtain 

V divergence = ~ 6{r = 0) JdrV- E(r) = -Nn <f(r = 0) , (3.56) 

where we used (3.37). This result for Vdivergence is merely a universal constant shift of 
energy, independent of the external potential V{r 2 ). It is proportional to N, whereas 
the leading behavior of energy eigenvalues is 0(N 2 ) ) as we will see momentarily. 
Thus, to leading order in j^, we shall neglect this term from now on. 

We should mention that this result for Vdivergence agrees, to leading order in 
with the one obtained in Eq. (11) in [23] for the quadratic potential V(r 2 ) = ^ r 2 . It 
was shown in [23] that this term was canceled against the divergent zero-point energy 
of the collective field 0(r). This cancellation should occur for any V(r 2 ), since our 
result (3.56) is manifestly independent of the external potential. 

Finally, the contribution of the external potential is 

V ext = N J dr<p(r)V(r 2 ). (3.57) 

Gathering all terms together, we obtain the collective hamiltonian as 

Hco11 = ^nI rfr Vn(r)-0(r)Vn(r)+^ J dr 0(r) E(r) 2 +N J dr 0(r) V(r 2 ) . (3.58) 

In order to display the large- N behavior of this theory explicitly, we rescale the 
coordinates as r = y/N^c. This is just the statement that TrM^M ~ iV / r 2 0(r)<ir is 
typically O(N), assuming <fi(\/Nx.) is supported in a disk |x| < O(N ). Consequently, 
0(r) = N^ 1 (p(x) (since 0(r) dr is invariant), and E(r) = y/NS(x) (from (3.35)), where 

£(*) = J dx'^f-^(x') (3-59) 
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is a quantity of O(N ). 

Also, from (2.14), V(r 2 ) = V(Ax 2 ) = iVV(x 2 ), where 

V(x 2 )=$>p(x 2 ) p . (3.60) 

p 

Finally, in order that the momentum dependent terms scale as the rest, we must 
rescale IT(r) = A 2 7r(x). The rescaled form of the hamiltonian (3.58) is thus 

Hcoii = iV 2 {|rfxQv7r(x)-v9(x)V7r(x) +^(x)^(x) 2 +^(x)V(x 2 )) 

- A(|dx^(x)-l)} , (3.61) 

where we have added a lagrange multiplier to impose the constraint (3.27). As ex- 
pected, the hamiltonian (3.61) and its eigenvalues, the ground state energy in partic- 
ular, scale as A 2 . 

The commutation relation of the rescaled canonically conjugate fields is 

b(x),vr(x')] = ^(5(x-x')-^) , (3.62) 

where we used (3.54). This commutation relation indicates that y?(x), 7r(x) and the 
dynamics of (3.61) become classical in the large- A limit. 

3.2.1 The Hamiltonian Equations of Motion: Fluid Dynamics Interpre- 
tation 

The Heisenberg equations of motion (as well as their classical analogs) resulting from 
(3.61) and (3.62) are easily found to be 

y9(x,t) = -iv-(Mx,t),V7r(x,t)}) , 

i.e., 

0(x, t) = -V • (y?(x, t) Vtt(x, f)) , (3.63) 



and 



*CM) = -i(Vvr(x,t)) 2 -j^fd^ [^(x)£(x) 2 + v,(x') (V(x' 2 ) - A) 



A 

+ L 2 
(3.64) 



31 



where 

A = JdK i(Vvr(x)) 2 + jdy J ^- ) (^(x) £(x) 2 + ^(x) (V(x 2 ) - A)) (3.65) 

is a constant. By assumption 13 , the coefficient in front of A in (3.64) tends to 
zero, and we shall neglect it. 

Eq. (3.63) just means that the constraint (3.27) is preserved by the dynamics. 

Let us take the gradient of (3.64) so as to make the equations of motion more 
symmetric. The 7r-dependent terms in the resulting equation can be lumped together 
into the combination <9iV7r+(<9 M 7r)V(<9 M 7r) = (9 t V7r + (V7r)-V(V"7r). Thus, the gradient 
of (3.64) may be written as 

d t X7ir + (Vtt) • V(Vtt) = -Vtc (3.66) 

where 

W (x; [tp]) = J dx! Q^(x') £(x') 2 + <p(x!) V(x' 2 ) - Av?(x')) . (3.67) 

Evidently, we should interpret the equations of motion as describing the isentropic 
flow of a two-dimensional eulerian fluid, with density </?(x, t), an irrotational velocity 
field v(x, t) = V7r(x, t), and enthalpy density u>(x; [</?]) given by (3.67). Eq.(3.63) is 
the continuity equation 

y?(x, t) + V • (y?(x, i)v(x, *)) = , (3.68) 

and (3.66) is just Newton's law 

<9 t v(x, t) + (v(x, t) ■ V) v(x, t) = -Vw (3.69) 

for isentropic flow, since in such flows the adiabatic variation of the enthalpy is related 

to the variation of the pressure P according to dw s = VdP = dP/ip. The fluid 

dynamical interpretation of the collective field equations of motion (3.63) and (3.64) is 

by no means special to our model. As we have already mentioned in the introduction, 

there is an analogous interpretation of the collective field of hermitean matrices in 

terms of a one-dimensional fluid [25, 26]. This is certainly a generic feature of field 

theoretic hamiltonians with a kinetic term of the form / rfx V7r(x) ■ </?(x)V7r(x) . 
13 See the discussion following (3.29). 
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4 The Ground State Energy and Distribution of 
Eigenvalues 



As was discussed in the previous section, <^(x, t), 7r(x, t) and the dynamics governed 
by (3.61) become classical in the large- N limit. More precisely, on general theoretical 
grounds [29], the quantum theory defined by (2.13) (and (3.61)) contains a special set 
of generalized coherent states, which may be used to construct a classical phase space 
and to derive a classical hamiltonian, whose resulting classical dynamics is equivalent 
to the large- N limit of the original quantum theory. The large- N ground state is then 
found by minimizing this classical hamiltonian. 

The solution of the classical equations of motion derived from (3.61) for (p(x.,t) 
and 7r(x, t), which minimizes H co u, should be interpreted as the expectation values 
of these fields in the ground state of (3.61). Evidently, the ground state (or vacuum) 
expectation value of 7r(x, t) should be null, and that of tp(x.,t) should be time inde- 
pendent. Thus, setting 7r (x, t) = (0|7r(x, t)\0) = and (0|</?(x, t)\0) = y?o( x ) into 
the classical equations of motion, we see that (3.63) is satisfied trivially, while (3.64) 
leads to the equation 

^/ rfx '(^(xW) 2 + ^x')V(x' 2 )-Ay,(x'))^ =0. (4.1) 

Eq. (4.1) is just the statement that y?o( x ) is the minimum of the effective potential 

Veff = JdK (^(x) £(x) 2 + y,(x) V(x 2 ) - A^x)) . (4.2) 

The external potential V(x 2 ) in our model is rotation-invariant. Thus, we would ex- 
pect that the eigenvalue distribution in the ground state should be rotation invariant 
as well. Furthermore, in a physically sensible model, the external potential V(x 2 ) 
must be confining, in order to balance the electric repulsion of the eigenvalues. The 
compromise between repulsion and confinement should produce a rotation-invariant 
charge density which is non-vanishing only inside a disk of finite radius R. 

In other words, it is sufficient to look for the minimum of (4.2) in the class of 
normalized density functions ip which depend only on r = |x|, and vanish outside a 
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disk of some radius R, which is fixed by the normalization condition 

Ft 



2vr 



/ <p(r) 
o 



r dr — 1 . 



(4.3) 



A radial charge density (p(r) is the source for a radial electrostatic field £(x) = £(r)r 
for which Gauss' law reads 



V-£=~(r£(r)) = 27r ¥ »(r), 



(4.4) 



and from which it follows that 



1 9 / 5£(r) 
r- 



2n5(r - r') 



(4.5) 



r dr \ 5ip(r') 

We stress that (4.5) holds only inside the support of <f(r), i.e., only in the region 
r < R. For r > R, the radial field 

5(r) = - (4.6) 
r 

identically, for any tp(r), rendering the variational derivative null in that region. 
Upon integrating (4.5) we obtain 

8£{r] 



6<p(r') 



2nr'9(r - r') 6(R - r) , 



(4.7) 



where we used the fact that £{r) arises only due to charges inside the disk of radius 
r to fix the integration constant. 

Thus, evaluating the effective potential (4.2) on the class of radially dependent 
densities, and taking the variation with respect to y?(r), we obtain, using (4.7), 

R 



sv f 



eff 



5<p(r) 



2irr 



\z 2 (r) + J 2it V {r')£{r') dr' + V(r 2 ) - A 



(4.8) 



Therefore, the desired minimum condition is 

R 



^£ 2 {r) + j 2n V {r')£{r') dr' + V(r 2 ) - A = . 



(4.9) 
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We can achieve considerable simplification of (4.9) by applying d r to both its sides 
and then use Gauss' law (4.4). We thus obtain the simple and elegant formula 



for the electric field inside the condensate of eigenvalues. 

By substituting the left-hand side of (4.4) for 2n(p(r') in (4.9), integrating, and 
using (4.10), we obtain an equation for the Lagrange multiplier A: 



Note from (4.10) that since S 2 (r) > 0, the eigenvalues can condense only in 
regions where V(r 2 ) is increasing (inside the disk r < R). This is obviously the result 
of electric repulsion among eigenvalues: in a region in which V(r 2 ) is decreasing, the 
external force pushes the eigenvalues away from the origin towrds larger radii, where 
they can keep farther apart from each other, until they reach a region in which V(r 2 ) is 
increasing and the external force is pointed inwards, which eventually confines them. 
This effect is demonstrated most clearly in the model in which the eigenvalues are 
confined inside a rigid disk of radius R, namely, V(r 2 ) = for r < R and infinite 
elsewhere. Evidently, in this case the eigenvalues are all stuck to the wall at r = R. 

A generic external potential V(r 2 ) might have several local wells inside the disk 
r < R, some (or all) of which might be occupied by eigenvalues. The eigenvalues 
will reside, of course, only in regions inside these wells in which V(r 2 ) is increasing. 
Let us assume that ^k(j') is an extremal eigenvalue density in which there are K 
occupied wells. Consider the kth. well (k = 1, . . . ,K). The eigenvalues will occupy 
in this well an annular region r]? m < r < r™ ax (and of course, = R). Thus, 

the eigenvalues in such a case will condense into a system of K concentric annuli, 
with voids in between the annuli. Inside the A;th annulus, the electric field is given by 
(4.10). In the void r™ ax < r < r™™ between the kth and k + 1st annuli, the electric 
field is obviously 



8 \r) = r-V(r 2 ) , 



r < R 



(4.10) 




(4.11) 




r 



k 



(4.12) 
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where v{r™ ax } = 2tt J (p(r)rdr is the fraction of eigenvalues confined up to r£ 



For k = K, (4.12) coincides with (4.6). It follows from (4.4) that £{r) is continuous 
wherever ip(r) is continuous. Thus, assuming V(r 2 ) has continuous derivatives, the 
electric field £(r) should be continuous throughout the plane. The K-annular ex- 
tremal eigenvalue density ip(r) depends upon the 2K parameters 14 r ™ n ' max , and the 
continuity conditions on the electric field at the 2K boundaries of the annuli provide 
the exact number of equations to determine these parameters. (The first of these 
conditions is of course £(r™ m ) = 0, and the last one E[r r ^ lx = R) = 1/R is equivalent 
to the overall normalization condition (4.3).) 

Given the external potential mentioned above, the extremum condition (4.9) might 
have several solutions v?R-(r), depending on the number K of occupied wells. In order 
to decide which of these extrema is the actual ground state configuration, we should 
evaluate the effective potential (4.2) at each of these extremal solutions and pick the 
the one with minimal energy. By tuning the couplings in V(r 2 ) we might induce quan- 
tum phase transitions between these possible multi-annular extremal configurations. 

In this paper we shall not pursue the generic case any further, with its multi- 
annular eigenvalue density configuration, and the fascinating possibility of inducing 
quantum phase transitions in this kind of matrix models. Instead, in what follows, we 
will solve for the ground state energy and eigenvalue distribution of the generic mono- 
tonically increasing external potentials, which necessarily yields a disk-like eigenvalue 
density configuration. However, in the explicit examples which will follow this expo- 
sition, we will compute the disk-annulus phase transition for the quartic potential. 

Thus, assume V(r 2 ) is such that y?o( r ) is voidless and supported in the disk r < R. 
The eigenvalues behave like positive point charges in our electrostatic desctiption of 
the matrix model. Thus, we should take the positive root of (4.10) for 



14 Note that once the 2K parameters r ™ m ' maa; are se t ; a ll the Vk are determined, since inside the 
filled annuli ip(r) is given by (4.4) and £{r) by (4.10). 



o 




r < R. 



(4.13) 
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From the continuity of £(r) at r = R and (4.6), we obtain 



RS(R) = 1, (4.14) 

which together with (4.13) forms an equation for R. ( Equivalent ly, we can derive this 
equation by substituting Gauss' law (4.4) into the normalization condition (4.3) .) 
Then, using (4.4), we can determine the eigenvalue density in the ground state as 

^^H^f- (4 - i5) 

From (4.11) and (4.14) we obtain 

A = ^ + V(i? 2 ) (4.16) 

for the Lagrange multiplier A. Finally, the ground state energy E is found by substi- 
tuting (4.13)-(4.15) (and 7r (x, t) = 0) in the collective hamiltonian (3.61), and using 
the normalization condition (4.3). After some algebra, we obtain 

e ° = N 'i dr (^ r£(r)) )(¥^ +v[r2) ) 

x ' 

4.1 Examples 

We end this paper by applying the formalism developed in this section in two explicit 
examples. 

4.1.1 The Quadratic Potential 

Consider the normal matrix model with quadratic potential (2.21), namely, 

V(r 2 ) = -m 2 r 2 (4.18) 
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with m 2 > 0. This is essentially the model studied recently in detail in [23], and our 
purpose here is simply to show that the general formalism developed above reproduces 
the results of [23] to leading order in 
From (4.13) - (4.15) we obtain 

S(r) = mr (4.19) 

and 

m 

<Po(r) = — = const. (4.20) 

7T 

in a disk of radius 

R = -$=. (4.21) 



'm 

Our result (4.20) is in agreement with Eq.(17) of [23] (with their A = 1). The ground 
state energy as found from (4.17) is 

E = ^N 2 m, (4.22) 
in agreement with (2.27) and also with Eq.(15) of [23], to leading order in 
4.1.2 The Quartic Potential 

Next we consider the quartic potential, which gives rise to a nonuniform ground state 
eigenvalue distribution y?o( r )- In this case the eigenvalues condense either into a disk 
or an annulus, depending on the sign of the coefficient of the quadratic term. The 
disk phase occurs for a positive coefficient of r 2 , and the annular phase occurs for a 
negative coefficient. 

(a) The Disk Phase 

Consider the quartic potential 

V(r 2 ) = imV + |r 4 (4.23) 
with m 2 ,g > 0. From (4.10), the electric field is 

£{r) = mr (l + ?£\ * (4.24) 
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and the density is 



Mr) = ~~U^, (4.25) 

in a disk of radius R, where R is obtained by solving (4.14). In the quartic case (4.14) 
leads to the cubic equation 

u 2 (u + l) = a 2 (4.26) 
for the positive dimensionless quantities 

u = ^ and a = At . (4.27) 

m 2 m 6 

The physical root of (4.26) should be chosen so as to have to correct value in the 
weak coupling and strong coupling limits. 

In the weak coupling limit a = ^ — > 0, the electric field (4.24) and eigenvalue 
density (4.25) tend manifestly to the corresponding quantities in the quadratic case. 
The root of (4.26) which tends smoothly to (4.21) is u ~ a « 1. 

The strong coupling limit a — > oo also yields simple expressions. The electric field 
tends to 

S(r) = ^r\ (4.28) 

the eigenvalue density tends to 

Mr) = ^r : (4.29) 

and the radius of the disk tends to 

R = g~*. (4.30) 

2 

The corresponding root of (4.26) is u ~ a» >> I. 

We shall not go into the solution of the cubic equation (4.26) for arbitrary coupling 
in any more detail. However, it is worth mentioning that the weak coupling and strong 
coupling branches of the solution correspond to two different roots of (4.26) each of 
which becomes real in its turn, which cross at a* = = 
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Finally, the ground state energy as found from (4.17) is 

mN 2 T9 _i ^ 9 ( , _5 ^ 1 

En 



Aa 2 



- ((« + 1)5 - l) - - ((« + 1)S - l) - - ((« + !)§-!) + ((„+ 1)* - l) 

(4.31) " 

It is easy to check that in the weak coupling limit (a,u « 1) this expression tends 
to (4.22), as it should. In the strong coupling limit (a,u » 1), it tends to 



E = 9 -^f(l + O(a~l)) , (4.32) 

which can be verified by plugging V(r 2 ) = |r 4 directly into (4.17). The scaling of the 
ground state energy in the strong coupling limit as is expected, of course, due to 
simple dimensional analysis. 

(b) The Annular Phase 

In the annular phase the quartic potential is 

V(r 2 ) = -i/iV + |r 4 , (4.33) 
with fx 2 ,g > 0. Then, from (4.10), the electric field is 

£{r) = ^r(^-iy (4.34) 

and the density is 

„2 



3 gr z _ i 

Mr) = ~ 1 ^—, (4.35) 



in an annulus r min < r < R. 

The radius of the inner boundary 

rLn = f (4-36) 

is found from the requirement that £ 2 (r) be positive. Note that near the inner 
boundary the density diverges as 

g 4 //A 2 1 



^ = — ij -^^ + 0(VF^-). (4.37) 

71 \ A J y / / min 
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The outer radius R is obtained by solving (4.14). The latter leads, in this case, to 
the cubic equation 

u 2 (u-l) = a 2 (4.38) 
for the positive dimensionless quantities 

u = ——r- and a = 4r , (4.39) 

in a similar manner to the disk phase. 

As in the disk phase, the physical root of (4.38) should be chosen so as to have 
the correct value in the weak coupling and strong coupling limits. 

In the annular phase, unlike the disk phase, the weak coupling limit a — — > 

is singular, since at g — 0, the potential (4.33) is unbounded from below and 
monotonically decreasing. The corresponding root of (4.38) in this limit is u — 

1 + a 2 + 0(a A ) which implies R = r min (l + + ^( fl4 ))- Thus, in the weak coupling 
limit the eigenvalues condense into a narrow annulus of very large radius ~ r m j„ = 

and width r min (i + 0(a 4 )) ~ \ (f)K 

The strong coupling limit a — > oo is essentially the same as in the disk phase, 

2 

as given by (4.28) - (4.30), and the corresponding root of (4.38) is u ~ as >> 1. 
In particular, note that in the strong coupling limit, the radius r min of the inner 
boundary shrinks like g~^a~^ = Ra~^, and the coefficient of the singularity (4.37) 
near the inner boundary diminishes like g^a~^, so that the voidless bounded density 
(4.29) will be obtained in the limit. 

We shall not go into the solution of the cubic equation (4.26) for arbitrary coupling 
in any more detail, except for mentioning the fact that (4.38) has only one real root, 
which interpolates smoothly between the weak and strong coupling limits. 

Finally, the ground state energy in the annular phase can be computed in a manner 
similar to the derivation of (4.17). One finds 
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which leads to the explicit expression 



Eq — 



UN 2 [9 




(4.41) 



4a 2 L7 



In the strong coupling limit (a,u >> 1), it tends to (4.32), namely the strong coupling 
limit of the disk phase. 

In the weak coupling limit (a << \,u ~ 1), on the other hand, (4.41) tends to 
the singular limit 



Thus, in the weak coupling limit, E tends to — oo, symptomatic of the fact that the 
model (4.33) does not exist at g — 0. 

(c) A Few Comments on the Disk-Annulus Phase Transition 

Starting from a given point (m 2 ,g) in the disk phase, with the potential (4.23), we 
can induce a quantum phase transition into the annular phase, by holding g fixed 
and tuning m 2 down all the way through m 2 = into negative values. The electric 
field £(r), eigenvalue density <fo(r) and ground state energy Eq are all continuous 
through the transition point at m 2 = jj? = 0. At the transition, they are given by the 
strong coupling formulas (4.28)-(4.30) and (4.32). Moreover, the annular eigenvalue 
distribution starts with a vanishing inner radius r min at the transition. 

This behavior is also expected of the disk-annulus phase transition in the zero 
dimensional (i.e., time independent) quartic ensemble of normal matrices. 

It is also interesting to compare the disk-annulus phase transition here with the 
corresponding transition in the zero dimensional quartic ensemble of complex matri- 
ces, studied in [8] and [9]. In the latter case, the Green's function (the resolvent) for 
the complex matrix M (i.e., the electric field, see (3.41)), the eigenvalue distribution 
in the complex plane, the Green's function for M^M and the density of singular values 
of M related to it, and thus the "free energy" of the Dyson gas of singular values, were 
all continuous through the transition. However, unlike the situation here, the disk 
phase there invaded well into the /z 2 > region in parameter space, and moreover, 




(4.42) 
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at the transition, the disk fragmented into an annulus of finite critical inner radius. 
(This occurred by having a central part of the disk, of finite size, progressively de- 
pleted, as the transition point was approached from the disk phse.) Also, the density 
of eigenvalues at the boundaries of either the disk or the annulus was always finite, 
unlike (4.37). 
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